Optical properties of graphene nanoribbons: the role of many-body effects 
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We investigate from first principles the optoelectronic properties of nanometer-sized armchair 
graphene nanoribbons (GNRs). We show that many-body effects are essential to correctly describe 
both energy gaps and optical response. As a signature of the confined geometry, we observe strongly 
bound excitons dominating the optical spectra, with a clear family dependent binding energy. Our 
results demonstrate that GNRs constitute ID nanostructures whose absorption and luminescence 
performance can be controlled by changing both family and edge termination. 
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Graphite-related nanoscale materials, such as 
fuUerenes and nanotubes, have long been the subject 
of an intense research for their remarkable proper- 
ties The recent discovery of stable, single-layer 
graphene 0, H, 3 has prompted the attention on a 
different graphitic quasi-lD nanostructure, i.e. graphene 
nanoribbons (GNRs). These systems have been the- 
oretically studied in the past decade s a a a i 

as 

simplified models of defective nanotubes and graphite 
nano- fragments. However, only very recently isolated 
nanometer-sized GNRs have been actually synthetized 
by etching larger graphene samples, or by CVD growth 
on suitably patterned surfaces [l^, EH, The pro- 

duction techniques advanced in these pioneering works 
are expected to become highly controllable, opening 
up new avenues for both fundamental nanoscience and 
nanotechnology applications. 

One of the most striking features of GNRs is the high 
sensitivity of their properties to the details of the atomic 
structure In particular, the edge 

shape dictates their classification in armchair (A), zigzag 
(Z) or chiral (C) ones, thus determining their energy band 
gaps. In addition to an overall decrease of energy gaps 
with increasing ribbon width, also observed experimen- 
. ^ ' tally 11| , theoretical studies predict a superimposed os- 
cillation feature [3, 14, [l^, which is maximized for 



A-GNRs. According to this behaviour, A-GNRs are fur- 
ther classified in three distinct families, i. e. N ~ 2>p~\, 
N = 3p, = 3p -|- 1, with p integer, where N indi- 
cates the number of dimer lines across the ribbon width. 
This fine sensitivity to the atomic configuration raise the 
opportunity to tailor the optoelectronic properties of A- 
GNRs by appropriately selecting both ribbon family and 
width. 

In spite of this interest, previous theoretical studies of 
the electronic (see e.g. Refs. [l^ 113) and optical prop- 
erties lJ| of GNRs were only based on the independent- 
particle approximation or on semi-empirical calculations. 
However, many body effects are expected to play a key 
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role in low dimensional systems 
due to enhanced electron-electron correlation 
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Moti- 



vated by this theoretical issue and by recent experimen- 
tal progress 13, U , l2| pursuing the potential of GNRs 
for nanotechnolgy applications, we have carried out ab 
initio calculations to study the effects of many-body in- 
teractions on the optical spectra of 1-nm-wide A-GNRs 
belonging to different families. 

In this Letter, we show that a sound and accurate 
description of the optoelectronic properties of A-GNRs 
must include many-body effects. We will demonstrate 
that there are many signatures of the non-local correla- 
tions occurring in these confined systems. First of all, 
quasiparticle corrections are found to be strongly state- 
dependent. Moreover, the optical response of A-GNRs 
is dominated by prominent excitonic peaks, with a com- 
plex bright-dark structure which would not have been 
even expected from an independent-particle framework. 
Both quasi-particle corrections and exciton binding ener- 
gies are found to exhibit an oscillating behaviour, accord- 
ing to the family classification. Finally, the electronic 
and optical properties of hydrogen passivated A-GNRs 
are compared with those of clean-edge ribbons: including 
many-body effects allows us to single out the impact of 
this edge modification on absorption and luminescence. 

The first-principles calculation of the optical excita- 
tions is carried out using a many-body perturbation the- 
ory approach, based on a three-step procedure [13] ■ As 
a preliminary step, we obtain the ground state electronic 
properties of the relaxed system, by performing a density- 
functional theory (DFT) supercell calculation, within the 
local density approximation (LDA) [i^, [ij]. Second, 
the quasiparticle corrections to the LDA eigenvalues are 
evaluated within the GqWo approximation for the self- 
energy operator, where the LDA wavefunctions are used 
as good approximations for the quasiparticle ones, and 
the screening is treated within the plasmon-pole approx- 
imation [ill. Third, the electron-hole interaction is in- 
cluded by solving the Bethe-Salpeter (BS) equation in 
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FIG. 1: (a) Optical absorption spectra of 1 nm wide 
hydrogen-passivated GNRs: N = 8 (1.05 nm wide), A'^ = 9 
(1.17 nm) and A'^ = 10 (1.29 nm). In each panel, the solid line 
represents the spectrum with electron-hole interaction, while 
the spectrum in the single-particle picture is in grey. All the 
spectra are computed introducing a lorentzian broadening, 
(b) Quasiparticle bandstructures. 



the basis set of quasielectron and quasihole states, where 
the static screening in the direct term is calculated within 
the random-phase approximation (RPA). Only the reso- 
nant part of the BS hamiltonian is taken into account 
throughout the calculations (Tamm-Dancoff approxima- 
tion), since we have verified that the inclusion of the 
coupling part does not affect significantly the absorp- 
tion spectra [1^. Moreover, only the case of light po- 
larized along the ribbon axis is examined, as a significant 
quenching of optical absorption is known to occur in ID 
systems for polarization perpendicular to the principal 
All the GW-BS calculations are performed with 
1, 



axis 



27| 



the code SELF 

To treat an isolated system in the supercell approach, 
we consider a separation of 40 a.u. between images in the 
directions perpendicular to the ribbon axis. Moreover, 
in both GW and BS calculations, we truncate the long- 
range screened Coulomb interaction between periodic im- 
ages, in order to avoid non-physical interactions 30]. Due 
to the rectangular geom etry of the system, we use a box- 
shaped truncation [31 1. 

We start by considering 1 nm wide hydrogen- 



passivated A-GNRs belonging to different families, 
namely = 8, 9, 10. Figure [1] (a) depicts their calcu- 
lated optical absorption spectra, while the quasiparticle 
bandstructures are shown in Fig. [1] (b) . All the results 
are summarized in Table HI The quasiparticle GW cor- 
rections open the LDA energy gaps at F by 0.72, 1.32 and 
1.66 eV for = 8, 9 and 10, respectively. These energy 
corrections are larger than those of bulk semiconductor 
with similar LDA gaps, due to the enhanced Coulomb 
interaction in low dimensional systems. In addition, a 
family modulation of the corrections can be noticed, with 
larger corrections for the GNRs with larger LDA gaps. 
The gap opening is accompanied by an overall stretch- 
ing of the banstructure of 17 — 22%, similar to the value 
found for graphene (about 20%) |32 |. 

In the absence of e-h interaction, such a bandstructure 
would result in the optical absorption spectra depicted 
in grey [Fig. [1] (a)], characterized by prominent 1-D van 
Hove singularities. The inclusion of the excitonic effects 
(solid black line) dramatically modifies both the peak po- 
sition and absorption line-shape, giving rise to individual 
excitonic states below the onset of the continuum, with 
binding energy of the order of the eV. 

The lowest-energy absorption peaks for A^ = 8 and 
9, labelled As and Ag, have the same character: in 
both cases, the principal contribution comes from optical 
transitions between the last valence and first conduction 
bands, localized in k-space near the F point [Fig. [T] (b)]. 
The binding energies for these lowest optically active ex- 
citons are 0.58 and 1.11 eV for A^ = 8 and 9, respectively. 
As compared to the first two systems, the A^ = 10 GNR 
shows a richer low-energy spectrum. Each noninteract- 
ing peak gives rise to a bright excitonic state [arrows 
and Bio in Fig- HI (b)], with binding energies of 1.31 and 
0.95 eV. In addition, the mixing of dipole forbidden tran- 
sitions between the same bands [arrows D in Fig[T](b)] is 
responsible for an optically inactive exciton degenerate in 
energy with Aiq . The D state thus provides a competing 
path for non radiative decay of optical excitations, which 
could affect the luminescence yield of the system. This 
feature results from transitions between pairs of bands 
very close in energy to each other, and is therefore ex- 
pected to be a common outcome for all A^ = 3p+l GNRs. 




FIG. 2: In-plane spatial distribution of the electron for a 
fixed hole position (black dot), corresponding to the lowest 
excitonic peak in the A'' = 9 case. The spatial density is 
averaged over the direction orthogonal to the ribbon plane. 
Dimension of the panel: 1.2 x 6.4 nm. 
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TABLE I: Energy gap (2nd and 3rd columns) and peak po- 
sition (4th column) for iV = 8, 9 and 10 A-GNRs, with (-H) 
and without hydrogen passivation of the edge sites. The rela- 
tive binding energies are reported in the last column. All the 
values are in eV. 



A further insight in the effects of electron-hole interac- 
tion is provided by the evaluation of the resulting spatial 
correlations. In Fig. [51 we plot the in-plane probabil- 
ity distribution of the electron for a fixed hole position 
(black dot), corresponding to the lowest excitonic state 
in the = 9 case. While the electron distribution ex- 
tends over the whole ribbon width, the modulation of 
the exciton wavefunction |V'(re;rh)|^ along the ribbon 
axis is entirely determined by the Coulomb interaction. 
Similar wavefunctions (not reported here) for the lowest 
excitons have been obtained for GNRs of different fami- 
lies, with spatial extentions [s^ of about 34, 23 and 18 
A for A^ = 8, 9 and 10, respectively. 

We now consider the case of clean-edge nanoribbons, 
since this simple variation of the structure has been of- 
ten suggested for ribbons obtained by high-temperature 
treatments or by dehydrogenation of hydrocarbons (sl-H. 
[s^l- This analysis allows us to further explore the role 
played by edge effects in the optoelectronic properties. 
Our results are summarized in Fig. [3] and Table [J As 
expected, the hydrogen removal leads to a major edge re- 
construction, with the appearence of carbyne-like struc- 
tures. In fact, the bond length for the edge dimers re- 
duces from 1.36 for the passivated ribbons to 1.23 A for 
the clean ones, pointing to the formation of C-C triple 
bonds at the edges. This edge modification leads to a 
variation of the energy gaps, such that the distinction 
between A^ = 3p — 1 and A'^ = 3£ families vanishes, in 
agreement with previous results [14| . 

In Fig. [3] (a), we report the quasiparticle bandstruc- 
ture for the A^ = 9 bare ribbon. The main difference 
with respect to its passivated counterpart is the pres- 
ence of edge-related bands (see arrows) in the low-energy 
optical window. Hence, we focus our attention on the 
properties of these edge states and their influence on the 
optical response. These states show the same energy dis- 
persion and real-space localization, irrespective of both 
family and size, already in the LDA framework 35|: due 
to this independence on bulk properties, their presence is 
reasonably expected for all non-passivated ribbons. The 
self-energy corrections to the LDA eigenvalues are similar 
to those of the passivated systems for the tt and it* bulk 
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FIG. 3: (a) Quasiparticle bandstructure of the A'' = 9 
hydrogen-free GNR. Arrows indicate the edge related single 
particle bands, (b) Plot of the GW quasiparticle energies 
vs the LDA energies, (c) Optical absorption spectrum, with 
(solid black) and without (grey) excitonic effects. The black 
arrow indicates the energy position of the optically forbidden 
edge-related exciton. Its excitonic wavefunction is depicted 
in panel (d), whose dimension is 1.0x2.2 nm. 



states. The edge states show quite a different correction, 
being deeper in energy and with a smoothed stretching 
with respect to the other bands [Fig[3] (b)]. This be- 
haviour is to be ascribed to the different degree of real- 
space localization between bulk and edge states, and it 
can be singled out by virtue of the non-local character of 
the self-energy operator in the GW framework, which is 
not correctly described within LDA. 

The aformentioned modification of the bandstructure 
results in a correspondent blueshift (A^ = 8) or redshift 
{N = 9) of the lowest excitonic peak, with Ag and Ag be- 
coming almost degenerate, with binding energies of about 
0.9 eV. For the case oiN = 10, we find an inversion of the 
first and second conduction bands, which results in the 
B^Q peak lying below Aj^q and D almost degenerate in 
energy with B^q. In addition, the edge states introduce 
an optically inactive exciton, which arises from transi- 
tions among several bulk valence bands and the conduc- 
tion edge states over the whole Brillouin zone. This edge 
exciton is present in all the studied nanoribbons and is 
located at about 1.4-1.7 eV (black arrow in Fig [3] (d)). 
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with very little dependence on family and size 35| . This 
results in the edge exciton being above the first excitonic 
peak for = 8 and 9, and between the first and the 
second peaks for N — 10. We remark that the accu- 
rate evaluation of quasi-particle corrections within GW, 
i.e. beyond the usual approximation based on a uniform 
band stretching on top of a rigid energy shift, is crucial 
to determine the exact energy position of the dark edge 
excitons relative to the bright ones. 

To better understand the character of the edge-related 
dark state, we plot its excitonic wavefunction for the case 
TV = 9 in Fig. [3] (d). The mixing of transitions over the 
whole Brillouin zone induces a strong localization of the 
edge exciton along the ribbon axis, with an extent of only 
'--^ 5 A, that is 4-7 times smaller than the Wannier-like 
bulk excitons (see Fig. [S]). 

In summary, we have found that the analysis of the 
electronic and optical features of GNRs requires a state- 
of-the-art approach within the many-body perturbation 
theory, and beyond the DFT framework. Many-body ef- 
fects reveal that nanosized A-GNRs retain a quasi-lD 
character, which induces the suppression of the van Hove 
singularity, typical of non-interacting ID systems, and 
the appearence of strong excitonic peaks in the optical 
absorption spectrum. The lowest excited states in GNRs 
are Wannier-like excitons and their binding energy as well 
as their luminescence properties are strongly dependent 
on the ribbon family. We investigate the role of many- 
body effects on the edge-states arising in non-passivated 
GNRs: our analysis could provide a practical tool for re- 
vealing the nature of the edges in realistic samples. We 
demonstrate that GNRs are intriguing systems with tun- 
able optoelectronic features, that we quantitatively eval- 
uate through our calculations. The present study calls for 
experiments addressing the optical response of GNRs: A 
combined theoretical and experimental understanding of 
ribbon size, family and edge-termination as control pa- 
rameters for their performance can be considered as the 
first step towards the design of graphene-based applica- 
tions in nanoscale optoelectronics. 
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